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We discuss the influence of large scale cosmic magnetic fields on the propagation of hadronic cosmic rays above 
10^^ eV based on large scale structure simulations. Our simulations suggest that rather substantial defiection 
up to several tens of degrees at 10^" eV are possible for nucleon primaries. Further, spectra and composition 
of cosmic rays from individual sources can depend on magnetic fields surrounding these sources in intrinsically 
unpredictable ways. This is true even if deflection from such individual sources is small. We conclude that the 
influence of large scale cosmic magnetic flelds on ultra-high energy cosmic ray propagation is currently hard to 
quantify. We discuss possible reasons for discrepant results of simulations by Dolag et al. which predict deflections 
of at most a few degrees for nucleons. We finally point out that even in these latter simulations a possible heavy 
component would in general suffer substantial deflection. 



1. Introduction 

The origin of ultra-high energy cosmic rays 
(UHECR) above lO^^eV (= lOEeV) has been 
challenging physicists for many years |1I2| . Sev- 
eral next-generation experiments, most notably 
the Pierre Auger experiment now under construc- 
tion and the EUSO project ^ are now trying 
to solve this mystery. 

Although statistically meaningful information 
about the UHECR energy spectrum and arrival 
direction distribution has been accumulated, no 
conclusive picture for the nature and distribu- 
tion of the sources emerges naturally from the 
data. There is on the one hand the approximate 
isotropic arrival direction distribution 5 which 
indicates that we are observing a large number of 
weak or distant sources. On the other hand, there 
are also indications which point more towards 
a small number of local and therefore bright 
sources, especially at the highest energies: First, 
the AGASA ground array claims statistically sig- 
nificant multi-plets of events from the same direc- 
tions within a few degrees |6I5| . although this is 



controversial and has not been seen so far by 
the fluorescence experiment HiRes The spec- 
trum of this clustered component is oc E~^'^ and 
thus much harder than the total spectrum |^. 
Second, nucleons above ~ 70EeV suffer heavy 
energy losses due to photo-pion production on 
the cosmic microwave background — the Greisen- 
Zatsepin-Kuzmin (GZK) effect [5] — which lim- 
its the distance to possible sources to less than 
~ 100 Mpc Pni- Heavy nuclei at these energies 
are photo-disintegrated in the cosmic microwave 
background within a few Mpc For a uniform 
source distribution this would predict a "GZK 
cutoff" , a drop in the spectrum. However, the 
existence of this "cutoff" is not established yet 
from the observations .12^. 

The picture is further complicated by the likely 
presence of large scale extra-galactic magnetic 
fields (EGMF) that will lead to deflection of any 
charged UHECR component. Magnetic fields are 
ubiquitous in the Universe, although their origin 
is still unclear l.'l . Magnetic fields in galaxies 
are observed with typical strengths of a few mi- 
cro Gauss. In addition there is some evidence 
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for fields correlated with larger structures such as 
galaxy clusters ^l]. Magnetic fields as strong as 
~ l/iG in sheets and filaments of the large scale 
galaxy distribution, such as in our Local Super- 
cluster, are compatible with existing upper limits 
on Faraday rotation |14I15| . It is also possible 
that fossil cocoons of former radio galaxies, so 
called radio ghosts, contribute significantly to the 
isotropization of UHECR arrival directions 

To get an impression of typical defiection an- 
gles one can characterize the EGMF by its r.m.s. 
strength B and a coherence length Ic- If we ne- 
glect energy loss processes for the moment, then 
the r.m.s. deflection angle over a distance r > Ic 
in such a field is e{E,r) ~ {2rlc/9y^^ /vl HT], 
where the Larmor radius of a particle of charge 
Ze and energy is ri, ~ E/{ZeB). In numbers 
this reads 

" {iW^ (lo^) ' 

for r > Ic- This expression makes it imme- 
diately obvious that fields of fractions of micro 
Gauss lead to strong defiection even at the high- 
est energies. This goes along with a time delay 
T{E,r) ~ re{E,df/4:, or 

r{E^r) 1.5xlOVrZ2(^^^) (2) 

^ (lOMpc) (ivV?) (icF^g) ' 

which can be millions of years. A source visi- 
ble in UHECR today could therefore be optically 
invisible since many models involving, for exam- 
ple, active galaxies as UHECR accelerators, pre- 
dict variability on shorter time scales. Strong 
deflection also limits the distance r a UHECR 
of given energy can travel during its energy loss 
time, which sometimes is called a "magnetic hori- 
zon" |18I19| . For example, if energy losses are 
included, for a space-filling field ~ 1 nC with 
= 1 Mpc the propagation distance of UHECR 
above lO^^eV is limited to - 400 Mpc For 
fields B ~ 0.1 ^G the magnetic horizon would be 
only ~ 20 Mpc 



2. Numerical Simulations: Large Scale 
Structure 

Quite a few simulations of the effect of extra- 
galactic magnetic fields (EGMF) on UHECR ex- 
ist in the literature, but usually idealizing as- 
sumptions concerning properties and distribu- 
tions of sources or EGMF or both are made: In 
Refs. I^OI^ll^^llAIS^ sources and EGMF follow 
a pancake profile mimicking the local supergalac- 
tic plane. In other studies EGMF have been ap- 
proximated in a number of fashions: as negligi- 
ble ^24,.2jT] . as stochastic with uniform statisti- 
cal properties |2HI27l28j . or as organized in spa- 
tial cells with a given coherence length and a 
strength depending as a power law on the local 
density Recently, Ref. |16I30| have carried 

out the first attempts to simulate UHECR prop- 
agation in a realistically structured universe; ad- 
ditional, independent calculations have now fol- 
lowed . So far all of these simulations have 
been limited to the case of nucleons. 

In Ref. the magnetized extragalactic en- 
vironment used for UHECR propagation is pro- 
duced by a simulation of the large scale structure 
of the Universe. The simulation was carried out 
within a computational box of 50 Mpc length 
on a side, with normalized Hubble constant h = 
i7o/(100 km s^i Mpc^i) = 0.67, and using a co- 
moving grid of 512^ zones and 256'^ dark matter 
particles. The EGMF was initialized to zero at 
simulation start and subsequently its seeds were 
generated at cosmic shocks through the Biermann 
battery mechanism Since cosmic shocks 

form primarily around collapsing structures in- 
cluding filaments, the above approach avoids gen- 
erating EGMF in cosmic voids. In this particular 
case (more are explored, see below) the resulting 
magnetic fields in collapsed structures are rather 
extended, as can be seen in Fig. ^ (top panel), 
significantly more than in the case of a uniform 
initial seed field (Fig. flower panel). The result- 
ing EGMF has been shown to be compatible with 
existing Faraday rotation measures with lines of 
sight both through clusters and the diffuse inter- 
galactic medium (I'H] . 

In Ref. [30] the following question was ad- 
dressed: which observer positions, and source dis- 
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Figure 1. Log-scale two-dimensional cuts through 
magnetic field total strength in Gauss (color scale 
in Gauss) for two EGMF scenarios. In the upper 
panel seed fields were injected at shocks, as in the 
scenario studied in Ref. |3U| , whereas in the lower 
panel, an initial uniform seed field was assumed. 
In both cases the seeds are chosen such that the 
final fields are of order ^ 1/i Gauss in a Coma-like 
cluster. For the latter case, fields in the voids were 
suppressed to minimize defiection. The observer 
is in the center of the figures and is marked by 
a star. The EGMF strength at the observer is 
~ 10~^^ G in the upper panel and ~ in the lower 
panel. The roughly spherically magnetized region 
about 3 Mpc to the lower right of the observer in 
the upper panel was used to simulate UHECR 
fluxes from an individual source, see Sect.|2| 



tributions and characteristics lead to UHECR dis- 
tributions whose spherical multi-poles for I < 10 
and auto-correlation at angles 9 < 20° are consis- 
tent with observations ? It was found that (i) the 
observed large scale UHECR isotropy requires the 
neighborhood within a few Mpc of the observer 
is characterized by weak magnetic fields below 
0.1 /iC, and (ii) once that choice is made, cur- 
rent data do not strongly discriminate between 
uniform and structured source distributions and 
between negligible and considerable deflection. 
Nevertheless, current data moderately favor a sce- 
nario in which (iii) UHECR sources have a density 
Hs ^ 10^^ Mpc^"^ and follow the matter distribu- 
tion and (iv) magnetic fields are relatively per- 
vasive within the large scale structure, including 
filaments, and with a strength of order of a fiG in 
galaxy clusters. A two-dimensional cut through 
the EGMF environment of the observer in a typ- 
ical such scenario is shown in Fig. ^ 

It was also studied in Ref. [301 how future data 
of considerably increased statistics can be used 
to learn more about EGMF and source charac- 
teristics. In particular, low auto-correlations at 
degree scales imply magnetized sources quite in- 
dependent of other source characteristics such as 
their density. The source characteristics can only 
be estimated from the auto-correlations halfway 
reliably if magnetic fields have negligible im- 
pact on propagation. Otherwise, the images of 
sources immersed in considerable magnetic fields 
are smeared out, which also smears out the auto- 
correlation function over several degrees. For a 
sufficiently high source density, individual images 
can thus overlap and sensitivity to source den- 
sity is consequently lost. The statistics expected 
from next generation experiments such as Pierre 
Auger ,3, and EUSO [l] should be sufficient to 
test source magnetization by the auto-correlation 
function [30! . 

In Ref. ^3L a constrained simulation which 
approximately reproduces the local universe on 
scales between a few Mpc and 115 Mpc from 
Earth was performed and the magnetic smoothed 
particle hydrodynamics technique was used to fol- 
low the EGMF evolution. The EGMF was seeded 
by a uniform seed field of maximal strength com- 
patible with observed rotation measures in galaxy 
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Figure 2. The cumulative distribution of deflec- 
tion angles a with respect to the line of sight 
to the sources of cosmic rays above 4 x 10^^ eV 
(upper panel) and above 10^" eV (lower panel). 
The blue (thick) and red (thin) curves showing 
larger and smaller deflections, respectively, are 
for the two EGMF scenarios from Fig. ^ upper 
and lower panel, respectively. The sources follow 
the baryon density and have average density — 
2.4 X 10~^ Mpc~'^. Shown are the averages over 
several realizations varying in the positions and 
luminosities Qi of individual sources, the latter 
assumed to be distributed as dris/dQi oc Q~^'^ 
with 1 < Qi < 100 in arbitrary units. 



clusters. 

The questions considered in this work were 
somewhat different. In Ref. [HJ deflections of 
UHECR above 4 x 10^^ eV were computed as a 
function of the direction to their source which 
were assumed to be at cosmological distances. 
Specific source distributions were not considered. 
The deflections typically were found to be smaller 
than a few degrees. 

Interestingly, however, there are considerable 
quantitative differences in the typical deflection 
angles predicted by the two EGMF models in 
Refs. |30l81j . These are not due to specific source 
distributions: In fact, for homogeneous source 
distributions, the average defiection angle for 
UHECR above 4 x lO^^ eV obtained in Ref. ^ is 
~ 61° above 4 x lO^^ eV and ~ 33° above 10^° eV, 
much larger than in Ref. |31j . 

In addition, even if the magnetic field strength 
is reduced by a factor 10 in the simulations in the 
environment of Fig. ^ upper panel, the average 
deflection angle is still ~ 28° above 4 x 10^^ eV, 
and ~ 10° above 10^*^ eV. This non-linear behav- 
ior of deflection with field normalization is mostly 
due to the strongly non-homogeneous character of 
the EGMF. 

Recently we have carried out a new simulation 
in which the initial magnetic seeds are provided 
by a uniform magnetic field (instead of the Bier- 
mann battery) to check how this different model 
affects our results. However, we find that typi- 
cal defiection angles change at most by a factor 
of 2. This is demonstrated by Fig. El which com- 
pares defiection angle distributions for the two 
EGMF scenarios shown in Figs. Q based on the 
Biermann battery model or uniform initial seeds, 
respectively. In both the above cases the mag- 
netic fields in a Coma-like cluster at 2; = were of 
order /iG. Also, in the case of uniform initial fields 
we have further suppressed fields in the voids to 
minimize deflections. 

Thus the residual differences in the predictions 
for the deflection angles of UHECR are proba- 
bly due to the different numerical models for the 
magnetic fields. Numerical issues may play an 
important role because they affect the amplifica- 
tion and the topological structure of the magnetic 
fields, both of which are important for the nor- 



5 



malization procedure, see below. A few examples 
are provided in the following. The resolution in 
Ref. is constant and in general better in fila- 
ments and voids but worse in the core of galaxy 
clusters than the (variable) resolution in Ref. |31| . 
Magnetohydrodynamic (MHD) codes are affected 
by numerical viscosity to a much larger extent 
than hydro codes; for this reason is it advanta- 
geous to evolve the magnetic field as a passive 
quantity because, while the amplification due to 
compression and stretching is fully accounted for, 
the additional dissipation introduced by the more 
complex MHD scheme is avoided. This is even 
more worrisome for the case of Ref. where 
the numerical dissipation is not quantified. It is 
also worth noting that the algorithm employed 
in Ref. does not guarantee a divergence-free 
magnetic field. While it is argued that in the high 
density regions (core of galaxy clusters) the varia- 
tions of the magnetic fields on a scale of a few res- 
olution elements are larger than the divergence- 
component, such conditions do not apply as easily 
in the low density, low resolution regions such as 
in filaments. Thus, one may wonder about the 
reliability of the numerical solution there. 

In conclusion, while in both simulations the 
magnetic fields are normalized to (or reproduce) 
the same "observed" values in the core of rich 
clusters, their values in the filaments are substan- 
tially different. While in Ref. [HOI the amplifica- 
tion in cluster cores may be underestimated, the 
magnetic fields in the filaments are not ruled out 
by available observational data ^1] . On the other 
hand, there is some concern for the reliability of 
the numerical results of Ref. PJj in the low den- 
sity regions, because numerical effects have not 
been quantified there. 

As a result, the magnetic fields obtained in 
Ref. jSni are considerably more extended than the 
ones in Ref. jSJ- This can be seen by comparing 
Fig. O with a similar plot in Ref. [33]: For the case 
shown in the upper panel of Fig. 13 about 10% of 
the volume is filled with fields stronger than 10 
nano Gauss, and a fraction of 10~^ is filled by 
fields above a micro Gauss. Corresponding filling 
factors for these field strengths in case of Ref. |i22| 
are < 10~* and < 3 x 10~^, respectively. This is 
the most likely reason for the much larger average 




Figure 3. The cumulative filling factors for 
EGMF strength in the two simulations shown in 
Fig. n above (decreasing curve) and below (in- 
creasing curve) a given threshold, as a function 
of that threshold. Upper and lower panels corre- 
spond to the panels in Fig. 
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deflections obtained in Ref. [201 as compared to 
the degree scale deflections discussed in Ref. . 

Since very little is currently known about the 
properties of intergalactic magnetic fields, the 
only solid conclusion that can be drawn at this 
stage is that the effects of EGMF on UHECR 
propagation are currently rather uncertain. 

Finally we note that these studies should be 
extended to include heavy nuclei 34. since there 
are indications that a fraction as large as 80% of 
iron nuclei may exist above 10^^ eV Prelim- 
inary results [Hi] give deflections of ^ 60° above 
10^*^ eV in our standard scenario shown in Figs.Q] 
and 121 upper panels. More importantly, even in 
the EGMF scenario of Ref. , deflections could 
be considerable and may not allow particle as- 
tronomy along many lines of sight: The distribu- 
tion of deflection angles in Ref. [HJ shows that 
deflections of protons above 4 x 10"'^^ eV of > 1° 
cover a considerable fraction of the sky. Suppres- 
sion of deflection along typical lines of sight by 
small filling factors of deflectors is thus unimpor- 
tant in this case. The deflection angle of any nu- 
cleus at a given energy passing through such ar- 
eas will therefore be roughly proportional to its 
charge as long as energy loss lengths are larger 
than a few tens of Mpc ■ Deflection angles of 
- 20° at ~ 4 X 10^^ eV should thus be the rule for 
iron nuclei. In contrast to the contribution of our 
Galaxy to deflection which can be of comparable 
size but may be corrected for within sufRciently 
detailed models of the galactic field, the extra- 
galactic contribution would be stochastic. Statis- 
tical methods are therefore likely to be necessary 
to learn about UHECR source distributions and 
characteristics. In addition, should a substan- 
tial heavy composition be experimentally con- 
firmed up to the highest energies, some sources 
would have to be surprisingly nearby, within a few 
Mpc, otherwise only low mass spallation products 
would survive propagation |37|. 

The putative clustered component of the 
UHECR flux whose fraction of the total flux 
seems to increase with energy may play a 
key role in this context. It could be caused by 
discrete sources in directions with small deflec- 
tion. Since, apart from energy losses, cosmic rays 
of same rigidity Z/A are deflected similarly by 



cosmic magnetic fields, one may expect that the 
composition of the clustered component may be- 
come heavier with increasing energy. Indeed, in 
Ref. it was speculated that the AGASA clus- 
ters may be consistent with consecutive He, Be- 
Mg, and Fe bumps. 

3. Numerical Simulations: Discrete 
Sources 

This brings us to the impact of EGMF on 
UHECR fluxes from discrete sources. In Ref. jSH] 
we investigated the impact on cosmic ray ob- 
servations above 10"'^^ eV of Mpc-scale magnetic 
fields of - 10"'^ G strength surrounding UHECR 
sources. The likely case was assumed that mag- 
netic fields within a few Mpc around Earth are in- 
significant. The trajectory simulations are based 
on the same large scale structure simulation as 
before, examples of whose EGMF distributions 
were shown in Figs. ^ and |21 upper panels. We 
find that such source fields can strongly mod- 
ify spectra and composition at Earth, especially 
for nearby sources for which the fields can con- 
siderably increase propagation times relative to 
both energy loss and photo-disintegration time 
scales and to the undeflected propagation time: 
Eq. (0) shows qualitatively that time delays at 
E ^ 10^" eV can easily reach ~ 10^ years for fields 
B - 10"'^ G extended over a few Mpc. This is in- 
deed larger than energy loss times ~ 3 x 10^ years 
and comparable to the straight line propagation 
time 10^ years. We found the following generic 
features: 

We note in passing that spallation and pion 
production could be enhanced around power- 
ful sources due to an increased infra-red back- 
ground . While we have not yet taken this into 
account in the present simulations, we expect this 
effect to be small compared to the EGMF effect. 

The spectra are considerably hardened relative 
to the injection spectrum at energies below the 
usual GZK-like cutoff where energy loss distance 
and source distance become comparable. This is 
caused by an interplay between diffusion and en- 
ergy loss: The flux of low energy particles is sup- 
pressed because diffusion spreads them out over 
a larger volume due to their much larger energy 
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Figure 4. Upper panel: Steady state all-particle 
spectra predicted by a scenario with a source at 
3.3 Mpc distance from Earth injecting iron pri- 
maries with a spectrum oc E^'^-'^ up to 4x 10^^ eV. 
The blue (thick) curves are for an EGMF of 
strength ~ 10~^ G surrounding the source and 
the red (thin) curve is without EGMF. Shown 
for comparison are the solid angle integrated 
AGASA ^ (dots) and HiRes-I 02] (stars) data. 
The solid straight line marks the injection spec- 
trum. All fluxes have been normalized at 60 EeV. 
Lower panel: The predicted composition in the 
same scenario. 



loss times. This is in contrast to the case of uni- 
formly distributed magnetic fields which in gen- 
eral lead to a steepening of the cosmic ray flux 
below the GZK cutoff. A hardened sub-GZK 
spectrum from individual sources would indeed 
be consistent with the hints of a hard clustered 
component in the AGASA data between 10^^ eV 
and 10^0 eV 0. 

Furthermore, for a nucleus of atomic mass A 
as injected primary, due to the kinematics of the 
photo-disintegration reactions a nucleon peak ap- 
pears at energy ^ E-^sx/A, where E-^^x is the 
maximal nucleus injection energy. This effect is 
the more prominent the harder the injection spec- 
trum. We also found that the details of spectra 
and composition depend significantly on the un- 
known details of the magnetic fields and the po- 
sition of the source therein and can thus not be 
predicted. 

An example demonstrating these effects is 
shown in Fig.^where a source of iron nuclei at 3.3 
Mpc distance to Earth is considered, once with- 
out any EGMF, and once with an EGMF concen- 
trated around the source and reaching ^ 10~^ G 
there. The observer is at the same position as 
indicated in Fig. ^ and the source is in the center 
of the roughly spherically magnetized region to 
the lower right of the observer in the upper panel 
of Fig. n The spectral modification for proton 
primaries (not shown) by the EGMF would be 
equally severe as in Fig. HI 

Further cutoffs towards low energies can be in- 
duced if the source is active only since a time 
smaller than the typical delay time at this en- 
ergy. Since the latter easily reaches ^ 10^ years 
at 10^^ eV in the scenario shown in Fig. 01 as can 
also be seen from Eq. jSJ, this is actually quite 
likely for active galactic nuclei which can have ac- 
tivity time scales below ~ 10* years. Obviously, 
the characteristics of such features also depend 
on unknown details of the source. 

Our simulations finally show that even for 
iron primaries, extra-galactic magnetic fields from 
large scale structure simulations are not strong 
and extended enough to explain the observed 
large scale isotropy of UHECR arrival direc- 
tions in terms of a single nearby source. This 
would require more homogeneous fields such as 
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Figure 5. The cumulative distribution of arrival 
direction off-sets from the source direction for 
UHECR above 4 x 10^^ eV in the scenario de- 
scribed in the text for iron primaries, correspond- 
ing to Fig. ^ (blue, thick curve) and proton pri- 
maries (red, thin curve). 



Next generation experiments such as the Pierre 
Auger Observatories jS] and the EUSO project 0] 
will accumulate sufficient statistics to establish 
spectra and distributions of composition and ar- 
rival directions from individual sources. A po- 
tentially strong influence of magnetic fields sur- 
rounding individual sources should thus be kept 
in mind when interpreting data from these exper- 
iments. This is true even if UHECR arrive within 
a few degrees from the source position. As Fig. [SI 
shows innij this is in fact the case for proton pri- 
maries in our example for a discrete source. 

4. Conclusions 

It thus seems evident that the influence of large 
scale cosmic magnetic flelds on ultra-high energy 
cosmic ray propagation is currently hard to quan- 
tify and may not allow to do "particle astronomy" 
along most lines of sight, especially if a signifl- 
cant heavy nucleus component is present above 




10^^ eV. However, even in our simulations, there 
are lines of sight along which deflection is only a 
few degrees, at least for proton primaries. Such 
directions might still be suitable for "particle as- 
tronomy" . 
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